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l/^ ■ We report the results of a phase-shift analysis (PSA) of the low-energy tt p elastic- 

^^N I scattering data. Following the method which we had set forth in our previous PSA 

[1], we first investigate the self-consistency of the low-energy Tr^p elastic-scattering 
databases, via two separate analyses of (first) the n^p and (subsequently) the -K~p 

(^ . elastic-scattering data. There are two main differences to our previous PSA: a) 

CN I we now perform only one test for the acceptance of each data set (based on its 

contribution to the overall x^) ^^d b) we adopt a more stringent acceptance criterion 
in the statistical tests. We show that it is possible to obtain self-consistent databases 

^ I after removing a very small amount of the data (4.57% of the initial database). 

H ' We subsequently fit the ETH model [38] to the truncated ir^p elastic-scattering 

" " " databases. The model-parameter values show reasonable stability when subjected 

to different criteria for the rejection of single data points and entire data sets. Our 
result for the pseudovector ttNN coupling constant is 0.0726±0.0014. We extract the 
scattering lengths and volumes, as well as the s- and p-wave hadronic phase shifts 
up to r = 100 MeV. Large differences in the s-wave part of the interaction can 
be seen when comparing our hadronic phase shifts with the current SAID solution 

~l/2 

(WI08); there is general agreement in the p waves, save for the 5i_ hadronic phase 

shift. 

PACS: 13.75.Gx; 25.80.Dj 
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1 Introduction 



This is tlie first of tliree papers dealing witli issues of tlie pion-nucleon (tiN) 
interaction at low energies (pion laboratory kinetic energy T < 100 MeV). 
The goal in this study is to update part of the material given in our previous 
phase-shift analysis (PSA) [1]. Hereafter, we will refer to that solution as 
UZH06, and to the one obtained in the present work as ZUAS12. We will 
show that the differences between these two solutions are small. In the second 
of the papers [2], we will address the self-consistency of the only ii^p elastic- 
scattering data which appeared in the meantime, i.e., of the measurements 
of Denz et al. (DENZ04) [3], give details on the problems we encountered 
in the analysis of these measurements, and provide evidence to support our 
decision to retain the UZH06 initial ir'^p elastic-scattering databases. In the 
third paper [4], we will analyse all available experimental data for the charge- 
exchange (CX) reaction n^p — > n^n and re-address the violation of the isospin 
invariance in the hadronic part of the vrA^ interaction. In our program, we 
make use of the electromagnetic (em) corrections of Refs. [5], which lead to 
the em- modified hadronic quantities [1]. 

We re-open the case of the ttN interaction at low energies because of three 
reasons. 

• There have been changes in the values of the physical constants which we 
had used in order to obtain our UZII06 solution. The largest of these changes 
concern the values of the charge radii of the pion and of the proton. All 
physical constants will be fixed here to the recommended values of the most 
recent report of the Particle-Data Group [6]. 

• In view of the fact that a large amount of n'^p elastic-scattering data became 
available in 2004, one might be under the impression that the important 
results and conclusions obtained with our UZH06 solution (which did not 
include these data) need to be revised; we will show that this is not the 
case. Due to numerous problems which we have encountered in the analysis 
of these data [2], we have decided not to include them in our database. 

• A large amount of CX data appeared after the completion of the UZII06 
PSA. In some cases, the corresponding reports of the experimental groups 
seem to cast doubt on the claims of isospin breaking [7,8]; therefore, it must 
be investigated whether the newly-obtained measurements invalidate those 
claims. 

We must stress that the physical quantities appearing here (i.e., the fit pa- 
rameters of Sections 3.1-3.4, the scattering lengths and volumes of Section 
3.4.2, the phase shifts of Section 3.4.3, etc.) are not purely-hadronic quanti- 
ties since they still contain residual em effects, which the stage-I em corrections 
of Refs. [5] do not remove. These residual effects relate in particular to the 



fact that, in principle, the (unknown) hadronic masses should be used (in the 
hadronic part of the vrA^ interaction) instead of the physical masses of the 
proton, of the neutron, and of the charged and neutral pion. Unfortunately, it 
is not possible at the present time to assess the significance of these effects. 
As a result, we must retain the cautious attitude of considering our physical 
quantities 'em-modified hadronic' (as we did in Ref. [1]). However, the repet- 
itive use of this term is clumsy. Therefore, we omit it, unless we consider its 
use necessary as, for instance, in Section 4 and in the captions of our tables 
and figures. 



2 Method 



2. 1 Formalism 



The determination of the observables from the hadronic phase shifts has been 
given in detail in Section 2 of Ref. [1]. For n^p scattering, one obtains the 
partial-wave amplitudes from Eq. (1) of that paper and determines the no- 
spin-flip and spin-flip amplitudes via Eqs. (2) and (3). Finally, the observables 
are evaluated from these amplitudes via Eqs. (13) and (14). For 7r~p elastic 
scattering, the observables are determined on the basis of Eqs. (15-20). 



2.2 Minimisation function and scale factors 



Similarly to our previous PSA [1], we make use of the minimisation function 
given by the Arndt-Roper formula [9]. The contribution of the j^^ data set to 
the overall x^ is: 
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where y^j^ denotes the i^^ data point of the j^^ data set, yjj the corresponding 
fltted value (also referred to as 'theoretical'), Sy^^^ the statistical uncertainty of 
the y^^^ data point, Zj a scale factor for the relative normalisation applying to 
the entire data set, Szj the corresponding uncertainty (reported or assigned), 
and Nj the number of data points in the data set. The fltted values yjj are 
obtained by means of parameterised forms of the s- and p-wave amplitudes. 
The values of the scale factor Zj are determined (for each data set separately) 
in such a way as to minimise x?- For each data set, there is a unique solution 



for Zj-. 

which leads to 
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The overall x^ = Y.jLi{Xj)min (where N stands for the number of data sets 
used in the fit) is a function of the parameters entering the modelling of the 
s- and p-wave amplitudes; these parameters were varied until x^ attained its 
minimal value Xmin- 

The part of {x'j)min which represents the pure random fluctuations in the 
measurements of the j*^ data set (i.e., the 'unexplained variation' in standard 
regression terminology) may be obtained from Eq. (3) in the limit 6zj — >■ oo, 
which is equivalent to removing the term {6zj)^'^ from the denominator of the 
second term on the right-hand side (rhs) of the expression; we will denote this 
value by (Xj)st- The variation which is contained in {x'j)min in excess of (x^)st 
must be associated with the contribution from the rescaling (floating) of the 
data set as a whole; the expression for {x'j)sc = {x'j)min — {Xj)st is 

, ^ i5zj)-'{T^,ySi:y^-yr)/iSy-ryy 



The scale factors which minimise only the first term on the rhs of Eq. (1) are 
obtained from Eq. (2) in the limit bzj — >• oo: 

. ^ T!:^.^y^li^y^f 



The scale factors of Eq. (2) are appropriate when investigating the goodness of 
the overall reproduction of a data set in terms of a reference solution (yielding 
the fitted values). On the contrary, those of Eq. (5) give maximal freedom to 
the baseline solution when determining the offset of a data set (with respect 
to that solution) and, as such, are more suitable when the emphasis is placed 
on the absolute normalisation, rather than on the overall reproduction, of a 
data set. 



The statistical uncertainty in the evaluation of the scale factors Zj is given by 

^^7 = , ^r ■ (6) 



When comparing the absolute normalisation of a data set with a reference so- 
lution, the statistical uncertainty of Eq. (6) must be taken into account, along 
with the normalisation uncertainty 6zj of that data set. When investigating 
the absolute normalisation of some specific data sets, namely of those with sus- 
picious absolute normalisation (e.g., see Section 3.4.6), the total uncertainty 
Azj = JSz'j + 5zj will be used. 



2.3 Statistical tests 



For a data set containing Nj data points, {Nj + 1) measurements have actually 
been made, the additional one relating to the absolute normalisation of the 
data set. Since the fit involves the fixing of each Zj at the value given in Eq. (2), 
the proper number of degrees of freedom (hereafter, the acronym DOF will 
stand for 'degree(s) of freedom', whereas NDF for the 'number of DOF'), 
associated with the j*^ data set, is just Nj. This implies that the quantity 
{Xj)min of Eq. (3) is expected to follow the x^ distribution with A'^- DOF. 

The essential difference between the present study and our previous PSA [1] 
is that only one statistical test for each data set will be performed here, the 
one involving its contribution {Xj)min to the overall x^'i iii R-^f. [1], we instead 
performed tests for the shape and for the normalisation of each data set. 

The p- value which is evaluated on the basis of {x'j)min and Nj will be compared 
with the confidence level Pmin for the acceptance of the null hypothesis (no 
statistically-significant effects); in case that the extracted p- value is below 
Pmin, the DOF with the largest contribution (to the {Xj)min value) will be 
eliminated in the subsequent fit. Only one point will be removed at each step, 
and the optimisation will be repeated. Data sets which do not give acceptable 
p- values (i.e., above Pmin) after the elimination of two of their data points 
(the absolute normalisation is also considered to be one data point) will be 
removed from the database. 

The second difference to Ref. [1], as far as the data analysis is concerned, 
relates to the choice of the confidence level which is assumed in the statistical 
tests; in Ref. [1], Pmin was set to about 2.70 • 10^^ (which is equivalent to a 
3a effect in the normal distribution). Herein, we will instead adopt the Pmm 
value which is associated with a 2.5a effect; this value is approximately equal 
to 1.24 • 10"^, that is, slightly larger than 1.00 • 10"^, which is 'commonly' 



(among statisticians) associated with the outset of statistical significance. In 
any case, only a few data (five DOF of the vr^p elastic-scattering databases) 
are affected by this more stringent acceptance criterion. 



3 Results 



The repetitive use of the full description of the databases is largely facili- 
tated if we adhere to the following notation: DB_|_ for the 7r~^p database; DB_ 
for the 7i~p elastic-scattering database; DB+/_ for the combined vr^p elastic- 
scattering databases. 

The initial DB^ comprises differential cross sections (DCSs) [10]- [18], analysing 
powers (APs) [19,20], partial-total cross sections (PTCSs) [21,22], and total 
(in fact total- nuclear) cross sections (TCSs) [23,24]. The initial DB_|_ consists 
of 364 data points, distributed among 54 data sets, 26 of which relate to the 
DCS, 3 to the AP, and 25 (all one- or two-point data sets) to the PTCS and 
TCS. The initial DB_ consists of 336 data points distributed among 36 data 
sets, i.e., 27 for the DCS ([12,13], [15]-[18], and [25]) and 9 for the AP ([19] 
and [26]- [28]). 

We now list the measurements which have not been included in our analysis. 

• The self-consistency of the n^p elastic-scattering DCS measurements of 
Ref . [3] will be addressed in detail elsewhere [2] . This is an enormous piece of 
experimental data, comparable in quantity to the database we established 
in our UZH06 PSA. Prior to their incorporation into a self-consistent set of 
data, the self-consistency of these measurements (as well as their compat- 
ibility with the established database) must be verified; in Ref. [2], we will 
come to a negative result. 

• The 70 7r~^p DCS measurements of Ref. [29] are obvious outliers in all ex- 
clusive analyses of the low-energy data; when using 'traditional' statistics, 
the inclusion of these values in a PSA is bound to introduce spurious effects 
(e.g., drifting of the parameters during the optimisation, entrapment of the 
minimisation algorithms in local minima, failing fits, etc.). We will further 
comment on these data in Section 3.4.6. 

• The 6 (3 for ir+p at 94.50 MeV, 3 for tt'p at 88.50 MeV) DCS measurements 
of Ref. [30] have not appeared in a form which would enable their straight- 
forward inclusion in our database. Furthermore, we are not convinced that 
the original DCS data may be retrieved by simply adding the contributions 
appearing in columns 4 and 5 of Table 1 of Ref. [30]. 

• Given that the 9 existing vr^p PTCSs and TCSs contain a component from 
CX scattering, they cannot be used; the inclusion of these data in any part 
of the analysis would perplex the discussion on the violation of the isospin 



invariance in the hadronic part of the ttN interaction [4] . 
The inclusion of the 28 AP measurements of Meier et al. (MEIER04) [31] in 
the fits would necessitate substantial modifications in the database structure 
and in the analysis software. This is due to two reasons: a) each of the three 
data sets, to which the measurements of Ref. [31] must be assigned, contains 
data taken at more than one beam energy and b) the last of the data 
sets contains measurements for both elastic-scattering reactions. Given the 
difficulty at present to include these measurements in our fits, we can only 
use them in testing the overall consistency of our approach across energies. 
(If a significant amount of 7r~p measurements had been acquired in that 
experiment, we could have investigated the consistency of our approach 
across reactions.) We will return to this subject in Section 3.4.5. 
We will also not use in the fits the scattering length obtained from the 
experimental result for the strong shift of the Is state in pionic hydrogen 
[32], after it has been corrected in Ref. [33] also taking into account the 
proper contributions of the 7n channel; the difference to the a'^'^ value of our 
UZH06 PSA (which is almost identical to the value we will obtain in this 
work) has been addressed in Ref. 



In order to give the data maximal freedom in the process of identifying the 
outliers, the two elastic-scattering reactions will be analysed separately using 
simple parameterisations of the s- and p-wave i^-matrix elements. The (small 
at low energies) d and / waves have been fixed herein from the current SAID 
solution (WI08) [34] . In the SAID analysis, the energy dependence of the d and 
/ phase shifts is determined from the region T > 100 MeV (i.e., from energies 
where these contributions are sizable). The largest of these phase shifts in the 
energy interval of this analysis, D15, does not exceed 0.27°. 

For the purpose of fitting, the standard MINUIT package [35] of the CERN 
library was used (FORTRAN version). Each optimisation was achieved on the 
basis of the (robust) SIMPLEX-MINIMIZE-MIGRAD-MINOS chain. All fits 
of the present work terminated successfully. 



3.1 Fits to the DB^ using the K -matrix parameterisations 



The parameterisation, which we will now describe, was introduced (and suc- 
cessfully applied to vr+p scattering) in Ref. [36]. For vr+p elastic scattering, the 
s-wave phase shift is parameterised as 

Qc cot 5o(^ = {aj^o+y^ + he + Cge^ , (7) 



where Qc and e are respectively the momentum and the pion kinetic energy in 
the centre-of-mass (CM) system. The pi/2-wave phase shift is parameterised 



according to the form 



Since the P3/2 wave contains the A(1232) resonance, a resonant piece in Breit- 
Wigner form is added to the background term, thus leading to the equation 

tan 5iV Qc = dsse + 6336 + ^ 2 ?i7i ' (9) 



with mA = 1232 MeV and Ta = 118 MeV [6]; qa is the value of the CM 
momentum at the resonance position. The third term on the rhs of Eq. (9) is 
the standard resonance contribution (e.g., see Ref. [37], p. 31). 

The use of the parametric forms in Eqs. (7)-(9), which do not impose any 
theoretical constraints (save for the expected low-energy behaviour of the K- 
matrix elements and for the Breit-Wigner form used to describe the A(1232)- 
resonance contribution), ensures that any outliers in the fits cannot be at- 
tributed to the inability of these forms to account for the energy dependence 
of the phase shifts; they are instead indicative of experimental problems. 

Following the procedure of Section 2.3, we found that, for the data sets of 
BRACK90 [16] at 66.80 MeV and JORAM95 [18] at 32.70 MeV, the p-values 
were still below Pmm after the removal of two data points (from each data set); 
as a result, these two data sets (with eleven and seven data points, respec- 
tively) were removed from the database. These two data sets have extremely 
low p-values and clearly stand out from the rest of the n^p data. Two data 
sets had to be freely floated, the BRACK86 [13] measurements at 66.80 and 
86.80 MeV. Finally, four additional single data points were removed. 

After the elimination of 24 DOF of the initial database, we obtain a truncated 
DB+ comprising 52 data sets with 340 DOF. The surviving data sets and the 
corresponding NDF are listed in Table 1; the numbers have been taken from 
the final fit to the truncated DB+ using the K-matrix parameterisations. 

Since seven parameters were used in order to generate the fitted values, the 
NDF in the fit to the initial database was 357; the minimal value of x^ w^-s 
675.7. For the truncated database with 333 DOF in the fit, the minimal value 
of x^ was 427.2. Therefore, the elimination of only 24 DOF of the initial 
database leads to an impressive decrease of the x^ by 248.5 units. At the same 
time, the p-value of the fit increased by over 17 orders of magnitude. 



3.2 Fits to the DB^ using the K -matrix parameterisations 



The / = 3/2 amplitudes, obtained in the final fit to the truncated DB_|_ using 
the i^- matrix parameterisations, were imported into the analysis of the DB_. 
In this part, another seven parameters were introduced, to parameterise the 
I = 1/2 amplitudes. As in the ii^p case, these parameters were varied pursuing 
the minimisation of the x^ function. Similar parametric forms were used as 
those given in Eqs. (7)-(9), with the parameters d^^. , 6i, ci, rfis, eis, du, and 

~l/2 

Cii. Of course, there is no resonance term in the expression for 6i_^_ ; instead, 

~l/2 

it is necessary to include the contribution of the Roper resonance in 6i_ : 

tan5}/Vge = due + e,,e' + ^^.^ /%,,, , (10) 



with 7727V = 1440 MeV and Fjv = 195 MeV [6]; q^ denotes the CM momentum 
at the Roper-resonance position. As we are dealing with energies below the 
pion-production threshold, Fjv is the elastic width. 

The DB_ was also subjected to the tests described in Section 2.3. As in the 
UZH06 PSA, the BRACK90 [16] 66.80 MeV data set (containing five data 
points) was marked as an outlier. Two additional single data points had to be 
removed. Furthermore, the WIEDNER89 [15] data set had to be freely fioated. 
After the elimination of 8 DOF of the initial DB_, we obtain a truncated 
database comprising 35 data sets with 328 DOF (see Table 2). Since seven 
parameters were used in order to generate the fitted values, the NDF in the 
fit to the initial database was 329; the minimal value of x^ was 528.1. For the 
truncated database with 321 DOF in the fit, the minimal value of x^ was 371.0. 
At the same time, the p-value of the fit increased by 9 orders of magnitude. 



3.3 Common fit to the DB^/^ using the K -matrix parameterisations 



Judged solely on the basis of the p- values, it appears that the truncated DB+ 
and DB_. are not of the same quality. However, there is an appropriate statis- 
tical measure for such a comparison of two quantities following the x^ distri- 
bution. In order to prove that the two databases are of different quality (that 
is, that they have not been sampled from the same distribution), the ratio 

F = 4^ (11) 



must be significantly different from 1. In this formula, the subscripts '+' and 
'— ' denote the two elastic-scattering reactions. The ratio F follows Fisher's {F) 



distribution. From the two final fits to the truncated databases using the K- 
matrix parameterisations, we obtain the score value of 1.110 for NDF^ = 333 
and NDF_ = 321 DOF in the two separate fits, which is translated into the 
p-value of 1.7 ■ 10^^. Therefore, the claim about the dissimilarity of the two 
databases cannot be sustained. As a result, it makes sense to analyse the two 
reactions in terms of a common optimisation scheme. 



In order to give the two elastic-scattering reactions equal weight, we multiplied 



{x^)min for each vr+p data set by 



w+ 



2^", 



and for each tt p elastic-scattering data set by 



W- 



2N_ 



where A^+ and AL represent the NDF in the two databases; we then added 
these quantities for all the data sets to obtain the overall x^ value. The appli- 
cation of these 'global' weights for the two reactions was made as a matter of 
principle; given that the values of A+ and N_ are very close, the effect of this 
weighting on our results is very small. 

The common fit to the truncated DB_|_/_ (detailed in Tables 1 and 2) was 
made, using the 14 parameters of the X-matrix parameterisations given in 
Sections 3.1 and 3.2. We did not find any additional data points (or data 
sets) which had to be removed; we concluded that the truncated DB_|_/_ are 
self-consistent. The common fit to the data yielded a x^ value of 792.4 for 654 
DOF. The set of the excluded DOF represents 4.57% of the initial database. 
In the following, we will use this truncated DB+/_. 



3.4 Common fits to the truncated DB^/_ using the ETH model 



The modelling of the hadronic part of the TriV interaction on the basis of 
the if-matrix parameterisations of Sections 3.1 and 3.2 is suitable as a first 
test of the self-consistency of the two elastic-scattering databases and as an 
efficient method for the identification of the outliers, yet neither does it provide 
insight into the underlying physical processes nor can it easily incorporate the 
important theoretical constraint of crossing symmetry. In order to accomplish 
these two tasks, we will next involve in the analysis a model based on Feynman 
diagrams, namely the ETH model. This model was introduced in Ref. [38] 
and was developed further throughout the early 1990s. The abihty of the 
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ETH model to account for the low-energy n^p elastic-scattering data has 
been convincingly demonstrated over the past two decades. 

The main diagrams on which this isospin-invariant model is based are graphs 
with scalar-isoscalar (/ = J = 0) and vector-isovector (/ = J = 1) t-channel 
exchanges, as well as the A^ and A s- and -u-channel graphs. The main contri- 
butions to the partial-wave amplitudes from these diagrams have been given in 
detail in Ref. [38]. The small contributions from the six well-established four- 
star s and p higher baryon resonances with masses up to 2 GeV have also been 
analytically included in the model; in fact, the only significant contributions 
come from the Roper resonance. The tensor component of the I = J = t- 
channel exchange was added in Ref. [8]; after this (insignificant) modification, 
no changes have been made to the model. 

The I = J = t-channel contribution to the amplitudes is approximated in 
the model by a broad tttt resonance, characterised by two parameters, G^r and 
Kcr- Its exact position has practically no effect on the description of the ttN 
scattering data or on the fitted values of G„ and K„; for a long time, it has 
been fixed at 860 MeV. The I = J = 1 t-channel contribution is described 
by the p-meson, with rrip = 775.49 MeV [6]; this contribution introduces two 
additional parameters: Gp and Kp. The contributions of the s- and w-channel 
graphs with an intermediate N involve the ttNN coupling constant q^nn and 
one additional parameter x representing the pseudoscalar admixture in the 
irNN vertex; for pure pseudovector coupling, x = 0. Finally, the contributions 
of the graphs with an intermediate A state introduce the coupling constant 
OttNA and the parameter Z (which is associated with the spin-1/2 admixture in 
the A field). The higher baryon resonances do not introduce any parameters. 

When a common fit of the ETH model to the data is made using all eight 
aforementioned parameters, it turns out that there is a strong correlation be- 
tween G^, Gp, and x; due to this correlation, it is not possible to determine 
the values of all three quantities. We have chosen to set x to 0; this choice 
is usually adopted in effective-field theoretical models of low-energy nN scat- 
tering. The common fits of the ETH model to the truncated DB+/_ will be 
performed on the basis of seven parameters: G^, K„, Gp, Kp, q-^nn^ dnNA, and 
Z. 



3.4-1 Model parameters 

The choice of the probability value below which data points must be excluded 
is difficult. We have adopted here the value of Pmin corresponding to a 2.5(T 
effect in the normal distribution. Recognising the subjectiveness in this choice, 
we consider that, in order to have confidence in the reliability of our analysis, 
it is necessary to verify that the fitted values of the seven model parameters 
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remain stable over a reasonably broad range of Pmm values; we followed the 
same strategy in Ref. [1]. Thus, in addition to Pmin ~ 1-24 • 10~^, the analysis 
was performed with a database reduced by using the Pmin values of about 
2.70 • 10~^ (equivalent to a 3a effect in the normal distribution) and 4.55 ■ 10~^ 
(equivalent to a 2cr effect in the normal distribution) lJ. 

Table 3 shows the values of the seven model parameters for the common fits to 
the truncated DB+/_ for the three selected values of Pmin- The uncertainties 
shown correspond to Pmm ~ 1-24 • 10~^. In fact, when the uncertainties are 



calculated with the Birge factor wx^/NDF included, they do not vary much 
with the value of Pmin- As p^m is increased, the truncated database fitted 



shrinks and so the raw uncertainties increase. However, the factor Wx^/NDF 
decreases as the fit quality improves (despite the decrease in NDF) and the 
two effects largely compensate. Table 3 shows that the results of the fit are 
reasonably stable as the criterion for the rejection of data points is varied L^. 

We see from Table 3 that the values ofo^ATAr are compatible; converted to the 
usual pseudovector coupling constant Lil, our result for Pmin ~ 1-24 ■ 10~^ is 

f^ = (J^] ^IeK = 0.0726(14) , 
47r \2mpJ Air ^ ^' 



where ^c denotes the mass of the charged pion and rup that of the proton. 
This result agrees well with the value we had obtained in our previous PSA 

The correlation (Hessian) matrix for the seven parameters of the ETH model 
is given in Table 4; the numbers correspond to the fit for Pmin ~ 1-24 • 10~^. 
This matrix, together with the uncertainties given in Table 3, enables the 
determination of predictions (and of their associated uncertainties) for the 
threshold constants, for the hadronic phase shifts and amplitudes, and for 
the observables at any combination of the energy and of the scattering angle. 



^ It must be mentioned that, in Ref. [1], the results for the different pmm levels 
had been obtained by applying cuts to the distribution of the residuals, as this 
distribution came out in the Pmin ~ 2.70 • 10"'^ solution. On the contrary, the 
analysis here is performed separately for the different Pmm levels. 
^ The effects, which are seen when increasing the Pmin value from the equivalent 
of a 2.5 to a 3a effect in the normal distribution, are due to the removal of the 
large JORAM95 vr+p data set at 44.60 MeV in the analysis of the data using the 
iC-matrix parameterisations of Section 3.1. A third data point from this data set 
(the measurement at 14.26°) must be excluded at the last step of the iteration, thus 
resulting in the removal of the entire data set; were this data set not excluded, there 
would have been almost no change in the parameter values shown in Table 3. 
^ Some authors redefine /^^^V' absorbing in it the factor Att. 
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Table 3 shows that the value of K^ is consistent with 0; the quality of the fit 
would deteriorate very little if this parameter were set to 0. The value of Go- 
is very little correlated with the values of the other five parameters. However, 
these parameters {Gp, Kp, q-^nni dnNA, and Z) are all strongly correlated with 
each other. (As expected, the correlations among the model parameters are 
smaller when no fioating of data sets is allowed in the fit.) 

Our results for the seven model parameters have shown remarkable stability 
over the years, from the period when the fits were performed to old, out- 
dated phase shifts (e.g., those of Refs. [39]) to the present times when the 
low-energy vrA^ measurements are directly fitted to. The database itself has 
changed significantly over the past twenty or so years, with the important 
contributions from experiments performed at the meson factories by differ- 
ent research groups, which made use of different apparatuses and techniques 
in their experiments. Our method of applying the em corrections has also 
changed. Finally, various approaches have been implemented in the optimi- 
sation (e.g., in the choice of the minimisation function, from 'standard' x^ 
functions [36], to robust statistics without any rejection of data [8], to the 
use of the Arndt-Roper formula [9] along with pruning the iiN databases [1]). 
The observed stability is indicative of the robustness one obtains in the results 
when involving the ETH model in the analysis of the low-energy ttN data. 



3-4-2 Threshold constants 

From the parameters of the ETH model and their uncertainties given in Table 
3 for Pmin ~ 1-24 • 10~^, as well as the correlation matrix given in Table 4, 
we determined the isoscalar and isovector s-wave scattering lengths and the 
isoscalar(isovector)-scalar(vector) p-wave scattering volumes. The results are: 

laV2 + ^a;;/.2 = 0.0033(12) /i,-i, 

'l 5o+ + I «o+ = -0.07698(60) /x'S 

I -at' + I at + \ -at + \ at' = 0-2039(19) /.-, (12) 

-i at + ^ ~at - \ ~at + \ -at = 0.1728(18) fi'^ 

^ at + I at - I at " ^ at = -0.1830(19) fif, 

'l at + I at + I at - ^ at = -0.06724(83) fif- 
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Converting these results to the famihar spin-isospin quantities, we obtain 



a, 



3/2 _ _nn7Q7n«^ ,,-1 ;;V2 

0+ - 



0.0737(16) fi~\ ool = 0.1573(13) fi'^ 



df^^ = 0.2090(20) /i;^ a}^^ = -0.03110(66) fi'^, (13) 

al^^ = -0.04124(78) /i;^ a}^^ = -0.0796(16) /i^l 

Our results for the s-wave scattering lengths ag^ and ao+ in Eqs. (13) are com- 
patible with those obtained in Refs. [36,8,1]; these values have been very stable 
over the last fifteen years. The large quantity of the elastic-scattering data be- 
low 100 MeV obtained at pion factories since 1980, when analysed without 
influences from the data obtained at higher energies, leads to results for the 
s-wave scattering lengths (and hadronic phase shifts) which are significantly 
different from those extracted via dispersion relations after also including the 
7r~p charge-exchange database in the analysis and using the measurements up 
to the few-GeV region. 

From the results in Eqs. (13), we obtain 
5^^ = I aj' + I ~al^^ = 0.0803(11) fi^' , 

in good agreement with the value extracted in Ref. [1]. We have already dis- 
cussed [1,33] the general disagreement (which, at present, is equivalent to 
a 4.6(7 effect in the normal distribution) of our d'^'^ value, extracted from 
Ti'^p elastic-scattering measurements, with the result obtained in the pionic- 
hydrogen experiments at the Paul Scherrer Institut (PSI) [32], after the cor- 
rections of Ref. [33] are applied. We comment further on this issue in Section 
6.1 of Ref. [4]. 



3.4-3 Hadronic phase shifts 

The s- and p-wave hadronic phase shifts, obtained from the common fit to the 
truncated DB_|_/_ using the ETH model, are given in Table 5. These phase 
shifts are also shown in Figs. 1-6, together with the current SAID solution 
(WI08) [34] and their five single-energy values (whenever available). 

It is evident from Figs. 1 and 2 that our values of the s-wave hadronic phase 
shifts (5o+ and 6q[^_ differ significantly from the SAID results. Our values of 
(5o+ are less negative, but converge towards the SAID values as the energy 
approaches 100 MeV; for 6q[^ , our values are consistently smaller. 

For the p-wave hadronic phase shifts 5i^ , 6i_ , and 6i\_ , inspection of Figs. 3-5 
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shows that there is general agreement between the two solutions (the differ- 
ences do not exceed about 0.1°). The significant difference in the p-wave part 
of the interaction occurs for 5iL ; our values are more negative. The values of 
5i_ , obtained in UZH06 and ZUAS12, are slightly different. This change is 
almost entirely due to the current use of a lower T^ value; at the time when 
the UZH06 PSA was carried out, the Particle-Data Group recommended the 
Pjv value of 227.5 MeV. 



Our phase-shift values are expected to be of interest in analyses involving the 
low-energy vrA^ interaction, as well as in the determination of the vriV a term 
(e.g., see Ref. [40]). 



We will now express our criticism concerning the SAID results at low energies. 



One has the impression that new ixN measurements enter the SAID database 
without regard of whether they comprise a self-consistent set and/or of 
whether they are at least marginally-compatible with the data which are 
already part of their database. Such a strategy could be less problematic, if 
they had implemented robust statistics in their data analysis; instead, their 
results are obtained with a 'standard' x^ function (i.e., the Arndt-Roper 
formula), and are thus expected to be sensitive to the presence of outliers 
in the database (in particular, to the presence of one-sided outliers). 
The SAID results at low energies are literally swamped by the measurements 
at higher energies. In case that the floating of the data sets is allowed (as is 
when using the Arndt-Roper formula in the optimisation), it is unavoidable 
that the low-energy behaviour of the ttN amplitudes will be influenced from 
higher energies. As a result, the low-energy experiments will be scaled sys- 
tematically in such a way as to match the trends of the amplitudes suitable 
for the higher-energy data. 

We have not found a published plot from the SAID group, showing the 
energy dependence of their scale factors below 100 MeV. 
The distribution of the normalised residuals and the function which is used 
in the optimisation are intimately connected. In case that a x^ minimisation 
function is used, the distribution of the normalised residuals must, in a self- 
consistent analysis, be a Gaussian centred at 0. Any effects observed in the 
distribution of the normalised residuals (e.g., significant offset, asymmetry 
of the distribution) are indicative of problems during the fitting procedure. 
We have not found this information in the SAID reports or in their web 
site. 

By forcing the data from all three reactions into an isospin-invariant analysis 
scheme, SAID cannot explore a possible violation of the isospin invariance 
in the hadronic part of the tiN interaction. 



15 



3.4-4 Scale factors and normalised residuals 

We will now comment on the distribution of the scale factors Zj obtained from 
the common fits using the ETH model. In a 'healthy' fit made on the basis 
of the Arndt-Roper formula, the data sets which must be scaled 'upwards' 
should (more or less) be balanced by those which must be scaled 'downwards'. 
Additionally, the energy dependence of the scale factors over the energy range 
of the analysis should not be significant. If these prerequisites are not fulfilled, 
the parametric forms used in the fits cannot adequately reproduce the data 
over the entire energy range. For both the vr+p (Fig. 7) and 7r~p (Fig. 8) 
elastic-scattering data sets, the values of Zj which lie above and below 1 roughly 
balance each other and there is no discernible energy dependence q|. Evidently, 
there is no subrange of the entire 30 to 100 MeV energy interval in which the 
data is better or worse fitted than for the rest of the range. 

A second issue which must be investigated in a fit involving the minimisation 
of any x^ function is whether the distribution of the normalised residuals Tj^, 
defined as 

^ ^^g- - y7 

(see Eq. (1)), evaluated on the basis of the optimal parameter values, is Gaus- 



nj 



g- 



sian Li|. In practice, one fits a Gaussian function to the distribution of the 
normalised residuals, i.e.. 



and investigates the quality of the fit (expressed through the corresponding 
X^ value and the NDF in the fit), as well as the asymmetry of the fitted 
distribution (e.g., expressed through the deviation of the extracted value of f 
from 0). The distribution of the normalised residuals is shown in Fig. 9, along 
with the optimal Gaussian function. The x^ value of this fit was 33.2 for 22 
DOF in the fit, whereas f = (1.9 ±4.7) • 10~^, i.e., compatible with 0. For the 
sake of completeness, we also give the optimal value and the uncertainty of 
the parameter B: B = 0.544 ± 0.038; the expectation value for B is 0.5. 



"^ The linear fit to the scale factors for the vr+p reaction yields the intercept of 
1.012 ± 0.018 and the slope of (-0.9 ± 2.4) • lO"'' MeV"^; the linear fit to the scale 
factors for the TT~p elastic-scattering reaction yields the intercept of 1.031 ± 0.018 
and the slope of (-4.5 ± 2.4) • 10"^ MeV-^. 

^ Of course, the values of {zj — l)/6zj, appearing in the second term on the rhs of 
Eq. (1), must also be included in the distribution of the normalised residuals. 
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3.4-5 Reproduction of the MEIERO4 measurements 

We will now discuss the MEIER04 measurements (which have not been in- 
cluded in our fits). We have created Monte-Carlo predictions for the AP corre- 
sponding to each of their 28 data points. For the three experimental data sets, 
the resulting values of Xmin were 12.5, 7.0, and 16.1, for 12, 6, and 10 DOF, 
respectively. The values of the scale factor for the three data sets (in the same 
order) are 1.011, 0.973, and 1.039; the reported normalisation uncertainty of 
the data is 3.5%. It is evident that our hadronic phase shifts reproduce the 
MEIER04 measurements very well; even the smallest p-value does not fall 
below 9.6 • 10^^. This is a good test of the consistency of our approach in 
describing the energy and angular dependence of the n^p elastic-scattering 
data. 



3.4-6 Reproduction of the BERTIN76 measurements 

We will now comment on the reproduction of the old measurements of Bertin 
et al. (BERTIN76) [29], which are 'traditionally' considered outliers in almost 
all modern PSAs. Reported in Ref. [29] were tt^p DCS measurements obtained 
in a broad angular interval, at seven beam energies (20.80, 30.50, 39.50, 51.50, 
67.40, 81.70, and 95.90 MeV). The experimental group did not report any nor- 
malisation uncertainties for their measurements; it cannot be excluded that, 
at those times, the absolute normalisation was not seriously investigated in 
the TriV experiments. 

We assigned rough normalisation uncertainties to the BERTIN76 data, on the 
basis of the values obtained from the modern low-energy experiments which 
properly reported this quantity, and analysed these measurements as if they 
comprised the entire DB_|_ at low energies. The analysis in terms of our general 
i^-matrix parameterisations of Section 3.1 revealed that the data set at 67.40 
MeV had to be eliminated due to its (very bad) shape. We were able to fit 
the remaining data successfully, and obtained the final x^ value of 55.0 for 53 
DOF in the fit. 

We subsequently investigated how well the phase-shift solution of the present 
study reproduces the BERTIN76 data; it failed. The basic problem with the 
BERTIN76 measurements lies with their absolute normalisation, not with 
their shape (though the shapes of the data sets at 39.50 and at 95.90 MeV 
do not pass our pmj„ criterion). For instance, for the data set at 20.80 MeV, 
Zj = 1.349 when the overall uncertainty Azj (defined at the end of Section 
2.2) is equal to 0.081; this is the most striking discrepancy in the data. The 
extracted Zj factors decrease almost linearly with T, reaching Zj = 1.027 at 
95.90 MeV. In view of these large energy-dependent effects in their absolute 
normalisation, we will continue excluding the BERTIN76 data in our PSAs. 
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4 Discussion 



In the present work, comprising tlie first of tliree papers dealing witli issues 
of the pion-nucleon (vrA^) interaction below pion laboratory kinetic energy of 
100 MeV, we report the results of a new phase-shift analysis (PSA) of the ii^p 
elastic-scattering data, using the electromagnetic (em) corrections of Refs. [5] . 

There are two main differences to the approach we set forth in our previous 
PSA [1], both pertaining to the method for the exclusion of outliers, single 
data points and entire data sets, in the optimisation phase. We now perform 
only one test for each data set on the only relevant quantity, namely on its 
contribution {Xj)min to the overall x^- In Ref. [1], we instead performed several 
tests (on shape, normalisation, etc.) for each data set; the use of only one test 
led to the exclusion of fewer data compared to Ref. [1] . The second difference 
concerns the imposition of a more stringent acceptance criterion of the null 
hypothesis in the statistical tests. Herein, we raised the minimal p- value (pmin) 
for the acceptance of the null hypothesis from the equivalent of a 3 to a 2.5a 
effect in the normal distribution; the latter value is closer to the 'common' 
choice (of most statisticians) as the outset of statistical significance. 

Similarly to Ref. [1], we first investigated the self-consistency of the low-energy 
ir^p elastic-scattering databases, via two separate analyses carried out (first) 
on the TT+p and (subsequently) on the 7r"p elastic-scattering data using sim- 
ple i^-matrix parameterisations. We found that it is possible to obtain self- 
consistent databases by removing the measurements of two n^p and one 7r~p 
data sets, as well as a few single data points; the removal of these outliers 
resulted in enormous reductions in the minimal x^ values for the separate fits 
to the two elastic-scattering databases using the i^-matrix parameterisations. 
We give all the details concerning the accepted data sets in Tables 1 and 2; 
these details may be useful in other analyses. The aforementioned results were 
obtained without imposing any theoretical constraints, save for the expected 
low-energy behaviour of the s- and p-wave i^-matrix elements and the Breit- 
Wigner form (see Eqs. (9) and (10)) for the contributions of the resonant 
terms. 

The ETH model of Ref. [38], based on s- and -u-channel diagrams with A^ 
and A in the intermediate states, and a and p t-channel exchanges, was sub- 
sequently fitted to the truncated combined Tr^p elastic-scattering databases. 
The model-parameter values showed reasonable stability when subjected to 
different criteria for the removal of data (see Table 3). Our result for the pseu- 
dovector ttNN coupling constant is 0.0726 ±0.0014. On the basis of the results 
of the model fits, we obtained the em-modified hadronic scattering lengths and 
volumes (see Section 3.4.2), as well as the s- and p-wave em-modified hadronic 
phase shifts up to T = 100 MeV (see Table 5 and Figs. 1-6). Large differences 



in the s-wave part of the interaction were found when comparing our hadronic 
phase shifts with the current SAID solution (WI08) [34] (see Figs. 1 and 2); 
there is general agreement in the p waves, save for the em-modified hadronic 
phase shift 6i_ . 

Apart from analysing our results in terms of the assumed confidence level in the 
statistical tests, we also investigated the possibility of a bias in the analysis. To 
this end, we examined the energy dependence of the scale factors Zj, shown in 
Figs. 7 (for the tt~^p data) and 8 (for the 7T~p elastic-scattering data), as well as 
the characteristics of the distribution of the normalised residuals (Fig. 9). We 
did not find any significant deviation for these quantities from the expectations 
in a successful optimisation. 
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Table 1 

The data sets comprising the truncated n^p database, the pion laboratory kinetic 
energy T (in MeV), the number of degrees of freedom (NDF)j for each data set, 
the scale factor Zj which minimises xi of Eq. (1), the values of {xi)min, and the 
p-value of the fit for each data set. The numbers of this table correspond to the 
final fit to the data using the i^-matrix parameterisations (see Section 3.1). 



Data set 


T 


(NDF), 


7 • 


\Xj jmin 


p-value 


Comments 


AULD79 


47.90 


11 


1.0146 


15.6800 


0.1534 




RITCHIE83 


65.00 


8 


1.0434 


17.2226 


0.0279 




RITCHIE83 


72.50 


10 


1.0047 


4.7383 


0.9080 




RITCHIE83 


80.00 


10 


1.0289 


19.0679 


0.0394 




RITCHIE83 


95.00 


10 


1.0327 


13.1452 


0.2157 




FRANK83 


29.40 


28 


1.0161 


19.4969 


0.8821 




FRANK83 


49.50 


28 


1.0458 


33.4861 


0.2183 




FRANK83 


69.60 


27 


0.9282 


23.6806 


0.6480 




FRANK83 


89.60 


27 


0.8614 


29.2304 


0.3498 




BRACK86 


66.80 


4 


0.8901 


2.4753 


0.6491 


freely floated 


BRACK86 


86.80 


8 


0.9380 


15.9483 


0.0431 


freely floated 


BRACK86 


91.70 


5 


0.9736 


11.6391 


0.0401 




BRACK86 


97.90 


5 


0.9723 


7.2220 


0.2046 




BRACK88 


66.80 


6 


0.9458 


11.5494 


0.0728 




BRACK88 


66.80 


6 


0.9554 


10.1689 


0.1177 




WIEDNER89 


54.30 


19 


0.9871 


14.7570 


0.7379 




BRACK90 


30.00 


5 


1.0830 


10.6037 


0.0598 


79.40° removed 


BRACK90 


45.00 


8 


1.0124 


7.6622 


0.4671 




BRACK95 


87.10 


8 


0.9730 


13.8367 


0.0861 




BRACK95 


98.10 


8 


0.9810 


14.8814 


0.0615 




JORAM95 


45.10 


9 


0.9600 


20.1704 


0.0169 


124.42° removed 


JORAM95 


68.60 


9 


1.0506 


8.2909 


0.5051 




JORAM95 


32.20 


20 


1.0138 


33.5992 


0.0290 




JORAM95 


44.60 


18 


0.9528 


29.9202 


0.0382 


30.74, 35.40° removed 


SEVIOR89 


98.00 


6 


1.0157 


5.4478 


0.4878 
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Table 1 continued 



Data set 


T 


(NDF), 


^j 


[Xjjmin 


p- value 


Comments 


WIESER96 


68.34 


3 


0.8899 


2.6924 


0.4415 




WIESER96 


68.34 


4 


0.9202 


3.8288 


0.4297 




KRISS97 


39.80 


1 


1.0129 


1.9961 


0.1577 




KRISS97 


40.50 


1 


1.0019 


0.1775 


0.6735 




KRISS97 


44.70 


1 


1.0027 


0.0768 


0.7817 




KRISS97 


45.30 


1 


1.0034 


0.0946 


0.7584 




KRISS97 


51.10 


1 


1.0246 


3.5130 


0.0609 




KRISS97 


51.70 


1 


1.0029 


0.0588 


0.8084 




KRISS97 


54.80 


1 


1.0077 


0.1738 


0.6768 




KRISS97 


59.30 


1 


1.0256 


1.2984 


0.2545 




KRISS97 


66.30 


2 


1.0497 


4.0101 


0.1347 




KRISS97 


66.80 


2 


1.0074 


0.5779 


0.7491 




KRISS97 


80.00 


1 


1.0136 


0.3366 


0.5618 




KRISS97 


89.30 


1 


1.0078 


0.2795 


0.5970 




KRISS97 


99.20 


1 


1.0560 


4.2590 


0.0390 




FRIEDMAN99 


45.00 


1 


1.0437 


2.3178 


0.1279 




FRIEDMAN99 


52.10 


1 


1.0182 


0.2772 


0.5986 




FRIEDMAN99 


63.10 


1 


1.0363 


0.4904 


0.4838 




FRIEDMAN99 


67.45 


2 


1.0517 


1.2397 


0.5380 




FRIEDMAN99 


71.50 


2 


1.0490 


0.8114 


0.6665 




FRIEDMAN99 


92.50 


2 


1.0429 


0.5872 


0.7456 




CARTER71 


71.60 




1.0921 


2.6734 


0.1020 




CARTER71 


97.40 




1.0498 


0.6952 


0.4044 




PEDRONI78 


72.50 




1.0121 


0.1329 


0.7155 




PEDRONI78 


84.80 




1.0311 


0.3258 


0.5682 




PEDRONI78 


95.10 




1.0230 


0.2030 


0.6523 




PEDRONI78 


96.90 




1.0167 


0.1330 


0.7153 
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Table 2 

The data sets comprising the truncated Tr~p elastic-scattering database, the pion 
laboratory kinetic energy T (in MeV), the number of degrees of freedom {NDF)j 
for each data set, the scale factor Zj which minimises Xj of Eq. (1), the values of 
{Xj)min, and the p-value of the fit for each data set. The numbers of this table 
correspond to the final fit to the data using the i^-matrix parameterisations (see 
Section 3.2). 



Data set 


T 


{NDF)j 


^i 


\Xj )min 


p-value 


Comments 


FRANK83 


29.40 


28 


0.9832 


30.9926 


0.3174 




FRANK83 


49.50 


28 


1.1007 


30.1075 


0.3581 




FRANK83 


69.60 


27 


1.0953 


25.6707 


0.5369 




FRANK83 


89.60 


27 


0.9479 


25.4255 


0.5506 




BRACK86 


66.80 


5 


0.9973 


13.9690 


0.0158 




BRACK86 


86.80 


5 


1.0036 


1.4172 


0.9224 




BRACK86 


91.70 


5 


0.9963 


2.8898 


0.7170 




BRACK86 


97.90 


5 


1.0002 


5.9408 


0.3120 




WIEDNER89 


54.30 


18 


1.1597 


23.6926 


0.1654 


15.55° removed, freely floated 


BRACK90 


30.00 


5 


1.0204 


4.8808 


0.4306 




BRACK90 


45.00 


9 


1.0536 


12.3031 


0.1968 




BRACK95 


87.50 


6 


0.9830 


10.3543 


0.1105 




BRACK95 


98.10 


7 


1.0092 


8.2223 


0.3134 


36.70° removed 


JORAM95 


32.70 


4 


0.9951 


4.0883 


0.3942 




JORAM95 


32.70 


2 


0.9527 


5.7974 


0.0551 




JORAM95 


45.10 


4 


0.9561 


12.4590 


0.0142 




JORAM95 


45.10 


3 


0.9462 


9.2581 


0.0260 




JORAM95 


68.60 


7 


1.0863 


14.2673 


0.0466 




JORAM95 


68.60 


3 


1.0314 


2.2747 


0.5174 




JORAM95 


32.20 


20 


1.0617 


21.3392 


0.3774 




JORAM95 


44.60 


20 


0.9462 


29.7408 


0.0742 




JANOUSCH97 


43.60 


1 


1.0420 


0.1682 


0.6817 




JANOUSCH97 


50.30 


1 


1.0364 


0.1557 


0.6931 




JANOUSCH97 


57.30 


1 


1.0830 


4.5370 


0.0332 
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Table 2 continued 



Data set 


T 


(NDF), 


^j 


{Xj)min 


p- value 


Comments 


JANOUSCH97 


64.50 


1 


0.9962 


0.0010 


0.9753 




JANOUSCH97 


72.00 


1 


1.3045 


4.8348 


0.0279 




ALDER83 


98.00 


6 


1.0335 


5.0919 


0.5321 




SEVIOR89 


98.00 


5 


0.9882 


1.5869 


0.9028 




HOFMAN98 


86.80 


11 


1.0020 


5.8362 


0.8841 




PATTERSON02 


57.00 


10 


0.9365 


11.0276 


0.3554 




PATTERS ON02 


66.90 


9 


0.9985 


4.3343 


0.8881 




PATTERS ON02 


66.90 


10 


0.9479 


18.0250 


0.0545 




PATTERS ON02 


87.20 


11 


0.9835 


8.2469 


0.6910 




PATTERSON02 


87.20 


11 


0.9945 


5.0493 


0.9288 




PATTERS ON02 


98.00 


12 


0.9954 


6.9731 


0.8594 





Table 3 

The values of the seven parameters of the ETH model obtained from the common 
fits to the truncated combined n^p elastic-scattering databases for three values of 
Pmin (the confidence level in the statistical tests); these three Pmm values correspond 
to a 3, 2.5, and 2a effect in the normal distribution, respectively. The uncertainties 
correspond to the fit for Pmin ~ 1-24 • 10~^. 





Pmin ~ 2.70 • 10-3 


p™„ « 1.24 • 10-2 


Pmin ~ 4.55 • 10-2 


uncertainty 


G^iGeV-^) 


27.43 


27.48 


27.37 


0.86 


K, 


0.014 


0.016 


0.075 


0.034 


G,{GeV-^) 


54.71 


54.67 


55.98 


0.61 


Kp 


0.66 


0.66 


1.35 


0.41 


gwNN 


12.84 


12.84 


13.07 


0.12 


9-kNA 


29.78 


29.77 


29.31 


0.26 


z 


-0.550 


-0.552 


-0.439 


0.056 
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Table 4 

The correlation matrix for the seven parameters of the ETH model for the common 
fit to the truncated combined n^p elastic-scattering databases for Pmin ~ 1.24-10^^. 





G. 


K. 


G, 


Kp 


gnNN 


g-TTNA 


Z 


G^ 


1.0000 


0.5095 


-0.0886 


-0.0378 


0.1030 


-0.1602 


-0.1999 


K. 


0.5095 


1.0000 


0.7314 


0.7974 


0.8847 


-0.9210 


0.7176 


Gp 


-0.0886 


0.7314 


1.0000 


0.9044 


0.9038 


-0.8487 


0.8977 


Kp 


-0.0378 


0.7974 


0.9044 


1.0000 


0.9522 


-0.9284 


0.9530 


g-TTNN 


0.1030 


0.8847 


0.9038 


0.9522 


1.0000 


-0.9497 


0.9216 


g-TTNA 


-0.1602 


-0.9210 


-0.8487 


-0.9284 


-0.9497 


1.0000 


-0.9001 


z 


-0.1999 


0.7176 


0.8977 


0.9530 


0.9216 


-0.9001 


1.0000 
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Table 5 

The values of the six s- and p-wave em- modified hadronic phase shifts (in degrees), 
obtained on the basis of the results of Tables 3 (for Pmm ~ 1-24 • 10^^) and 4. 



T{MeV) 


6li' (S31) 


6li' (Sll) 


~5li' (P33) 


6f_^ (P31) 


~6li' (P13) 


~6lL' (Pll) 


20 


-2.375(34) 


4.189(27) 


1.2787(94) 


-0.2239(45) 


-0.1588(37) 


-0.3687(80) 


25 


-2.772(36) 


4.673(29) 


1.817(12) 


-0.3083(63) 


-0.2153(52) 


-0.486(11) 


30 


-3.164(37) 


5.105(30) 


2.431(15) 


-0.3996(83) 


-0.2747(68) 


-0.602(14) 


35 


-3.555(37) 


5.496(31) 


3.122(18) 


-0.497(11) 


-0.3361(86) 


-0.714(17) 


40 


-3.949(37) 


5.852(33) 


3.892(21) 


-0.599(13) 


-0.399(11) 


-0.820(20) 


45 


-4.345(37) 


6.180(34) 


4.744(23) 


-0.706(16) 


-0.463(13) 


-0.918(23) 


50 


-4.746(37) 


6.482(37) 


5.683(24) 


-0.816(18) 


-0.527(15) 


-1.006(27) 


55 


-5.151(38) 


6.760(39) 


6.715(26) 


-0.931(21) 


-0.591(17) 


-1.084(30) 


60 


-5.561(39) 


7.018(42) 


7.845(27) 


-1.049(25) 


-0.655(20) 


-1.150(34) 


65 


-5.977(41) 


7.256(46) 


9.081(29) 


-1.169(28) 


-0.719(23) 


-1.204(38) 


70 


-6.397(44) 


7.476(51) 


10.433(31) 


-1.293(32) 


-0.782(26) 


-1.244(42) 


75 


-6.823(48) 


7.679(56) 


11.909(35) 


-1.419(36) 


-0.844(29) 


-1.271(46) 


80 


-7.254(53) 


7.865(61) 


13.519(40) 


-1.547(41) 


-0.906(32) 


-1.283(50) 


85 


-7.690(59) 


8.036(67) 


15.277(48) 


-1.678(45) 


-0.966(36) 


-1.281(55) 


90 


-8.131(67) 


8.192(74) 


17.193(59) 


-1.811(50) 


-1.026(39) 


-1.263(60) 


95 


-8.577(75) 


8.334(81) 


19.282(73) 


-1.946(55) 


-1.084(43) 


-1.229(66) 


100 


-9.028(85) 


8.462(89) 


21.556(90) 


-2.083(61) 


-1.141(48) 


-1.178(71) 
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Fig. 1. The em-modified hadronic phase shift 6q_^ (S31) from the present work (sohd 
curve); the band around our solution indicates Icr uncertainties. Shown also is the 
current SAID solution (WI08) [34] (dashed curve), along with their five single-energy 
values (at T = 20, 30, 47, 66, and 90 MeV). 
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Fig. 2. The em-modified hadronic phase shift (5q^ (Sll) from the present work (sohd 
curve); the band around our solution indicates Icr uncertainties. Shown also is the 
current SAID solution (WI08) [34] (dashed curve), along with their five single-energy 
values (at T = 20, 30, 47, 66, and 90 MeV). 
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Fig. 3. The em-modified hadronic phase shift 6i_^ (P33) from the present work 
(solid curve); the band around our solution indicates la uncertainties. Shown also 
is the current SAID solution (WI08) [34] (dashed curve), along with their five sin- 
gle-energy values (at T = 20, 30, 47, 66, and 90 MeV). To enable the meaningful 
comparison of the values contained in this figure, an energy-dependent baseline 6r 
(= (0.20 ■ T + 1.54)r • 10-2, with T in MeV and 6r{T) in degrees) was subtracted 
from all data. 
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Fig. 4. The em-modified hadronic phase shift (5^_ (P31) from the present work (soHd 
curve); the band around our solution indicates la uncertainties. Shown also is the 
current SAID solution (WI08) [34] (dashed curve). 
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Fig. 5. The em-modified hadronic phase shift 5^^ (P13) from the present work (solid 
curve); the band around our solution indicates \a uncertainties. Shown also is the 
current SAID solution (WI08) [34] (dashed curve). 
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Fig. 6. The em-modified hadronic phase shift (5^_ (Pll) from the present work (soHd 
curve); the band around our solution indicates la uncertainties. Shown also is the 
current SAID solution ( WI08) [34] (dashed curve) , along with their five single-energy 
values (at T = 20, 30, 47, 66, and 90 MeV). 
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Fig. 7. The scale factors Zj for the Tr~^p data, obtained from the common fit to 
the truncated combined n^p elastic-scattering databases using the ETH model (see 
Section 3.4). The values, corresponding to the two data sets which were freely floated 
(see Table 1), have not been included; in the case of free floating, (x?)sc = 0. 
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Fig. 8. The scale factors Zj for the ■7r~p elastic-scattering data, obtained from the 
common fit to the truncated combined n^p elastic-scattering databases using the 
ETH model (see Section 3.4). The value, corresponding to the data set which was 
freely floated (see Table 2), has not been included; in the case of free floating, 
{X%c = 0. 
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Fig. 9. The distribution of the normahsed residuals, obtained from the common fit 
to the truncated combined vr^p elastic-scattering databases using the ETH model 
(see Section 3.4). Also shown (solid curve) is the optimal Gaussian fit to the data. 
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